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The excited state energy spectrum of a two-dimensional dimer system SrCu2 (603)2 has been 
investigated by high field ESR. Various types of multiple-triplet bound states-excited singlet, 
triplet and quintet states are identified, besides two non-degenerated one-triplet excitations. 
Extremely strong intensity is found for these higher order bound states, which is an important 
consequence of the strong suppression of the one-triplet process due to the orthogonality of 
dimers. A new bound singlet state is found slightly below the energy gap of the one-triplet state, 
indicating the proximity to a quantum critical point. A band of two-triplet bound state splits 
into at least eight sharp and discrete levels, which indicates the contribution of distant neighbor 
interactions. A bound quintet state is observed, which is the direct evidence of two-triplet object 
made up of four 5 = 1/2 spins. The energy of the lowest quintet becomes lower than that of the 
one-triplet excitation before the closing of the lowest spin gap. It indicates that the two-triplet 
bound state touches to the ground singlet state before the one-triplet state does. The present 
investigation exhibits an essential role of various multiple-triplet bound states in the magnetic 
excitation spectrum of SrCu2 (603)2- 

KEYWORDS: SrCu2 (603)2, spin gap, ESR, magnetic excitation, high magnetic field 



§1. Introduction 

Coupled chain and dimer systems have been attracted 
much attention in this decade. One of the interests 
is in a complex and fruitful crossover from a well- 
known purely one-dimensional antiferromagnet to a two- 
dimensional(2D) antiferromagnet, i.e., the basis of a 
high-temperature superconductor. Experimentally, a 
2D spin gap system is quite rare, in contrast with a rapid 
growth of new quasi-one dimensional compounds. One 
of fundamental difficulties in studying the 2D system ex- 
perimentally is that it tends to exhibits three dimen- 
sional(3D) long range order. However, there is an effec- 
tive way to prevent the 3D-order:introduction of frustra- 
tion into a system. One of the most successful exam- 
ples is SrCu2(B03)2:a 2D dimer system with consider- 
ably strong frustration.'^) 

In this compound, a layer of S* = 1/2 Cu^+ ions is 
sandwiched by layers of Sr ions in a tetragonal unit cell.^-* 
A dimer unit is composed of a neighboring pair of edge- 
shared planar rectangular CUO4 and these dimers con- 
nect orthogonally by way of a triangular planar BO3, 
providing a unique 2D-network as shown in Fig. 1. This 
2D lattice is topologically equivalent to a 2D square lat- 
tice with additional alternating diagonal interactions, for 
which the direct product of the singlet pairs is the exact 
ground state of the system, as proven by Shastry and 
Sutherland.**) 

Although the ground state is known rigorously as a 
dimer solid, the excited state is expected to be very rich. 
A low- lying excited triplet, which is dominant for spin 



gap compounds in general, is considerably suppressed 
due to the extreme localization. Higher order multiple- 
triplet states are instead stabilized by the kinetic energy 
called "correlated hopping". Those interesting features 
of the excited state have been studied experimentally and 
theoretically, although a complete understanding has not 
been obtained so far.^"^^) 

Another important issue is that this compound is lo- 
cated very close to a quantum critical point. The ratio 
J2/J1 is evaluated to be 0.63-0.68, where Ji and J2 are, 
respectively, the intradimer and interdimer interactions 
as shown in Fig. 1 (a)."'^'^-' The value is just below the crit- 
ical value of J2I Ji=Q-7G between the dimer state and the 
Neel-ordered state. More recently, it has been proposed 
that a novel ground state, a resonating plaquette singlet 
state, exists in between these two states.**-' Although 
the present compound is in the dimer phase, a precursor 
effect of other types of ground states may appear, espe- 
cially in excited states, for the proximity to the quantum 
critical point. 

An interesting feature also appears in the magne- 
tization curve with quantized plateaux at one-third, 
one-quarter and one-eighth of the Cu saturated mo- 
ment.*^' These plateaux are caused by the extreme 
localization of the triplet excitation. Namely, the excited 
triplets cannot propagate freely in the 2D lattice owing 
to the orthogonal dimer arrangement and consequently 
organize a regular lattice composed of spatially separated 
triplets when the density of the triplets is commensurate 
with the underlying lattice. 

As discussed above, it is very interesting to study the 
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excited state spectrum of the present system, especially 
as a function of a magnetic field. Among many differ- 
ent methods, ESR has established its unique status as a 
probe of magnetic excitations, particularly in studying 
the high field property.^^"^'*) In fact, the one-triplet as 
well as the multiple-triplet excitations were reported in 
our previous paper. However, the richness of the spec- 
trum and the limited frequency range prevented us to 
complete the experimental determination of a full spec- 
trum. In the present paper, a nearly complete excited 
state energy spectrum of the SrCu2 (603)2 up to the two- 
triplet continuum and up to 40 T is obtained. 

The format of the paper is as follows. After briefly 
mentioning the experimental procedure, we show details 
of experimental results. The energy diagram, selection 
rule and the symmetry of one-triplet and multiple-triplet 
excitations are examined. The anomalous behavior of 
the magnetic excitation and its relation to the magneti- 
zation around the critical field is discussed in terms of 
the localization of the bound states. 



Fig. 1. 2D-dimer plane viewed along the c-axis with respect 
to z=0, or 0.5. Open, shadowed and closed circles, denotes 
Cu, O and B ions, respectively. The size of a circle indi- 
cates the shift along the c-axis(bigger/smaller, corresponding 
upward/downward). An arrow on Cu shows a tilting of the 
principle axis of g-tensor. A and B distinguish two types of 
dimers(with dashed ellipsoids) which are connected orthogonally 
each other. Thin square shows a unit cell. Ji (thick dashed line) 
and J2(thin dashed line) are the nearest neighbor and the next- 
nearest-ncighbor interaction in the plane, respectively. The inset 
is the [l-10]-axis view of the A-dimer. Ellipsoids schematically 
represent the anisotropy of g-tensor. 



§2. Experimental procedure 

Submillimctcr wave ESR measurements have been per- 
formed up to 1.3 THz and in pulsed magnetic fields up 
to 40 T. An optical pumped far-infrared laser, backward- 
travelling wave tubes and Gunn oscillators have been em- 
ployed as radiation sources. We have employed a simple 



transmission method with Faraday configuration where 

the propagation vector of the incident radiation is aligned 
parallel to the external magnetic field. The polarization 
of light is random in the plane normal to the propagation 
vector. An InSb is used as the detector. The detail of 
the ESR system is given in the references. ^^'^^^ 

High-purity bulk single crystals of SrCu2(B03)2 were 
grown by the travelling solvent floating zone (TSFZ) 
method from a polycrystalline SrCu2 (603)2 and a sol- 
vent LiB02 using FZ-TIOOOON lOKW high pressure type 
(Crystal System Inc.) with four halogen lamps as heat 
sources. Oxygen gas (Po2=l atom, 99.999 %) was being 
flowed during the growth process. By means of Lauc 
X-ray back-reflection, the grown materials were checked 
and the crystal axes were determined. A piece of the sin- 
gle crystal with the dimensions of about 4 mmx4 mmx 1 
mm was used for the ESR experiments. 

§3. One-triplet excitation 

In the present work, we mainly focus on ESR spectra 
at 1.6 K and below the critical field He at which the 
gap closes. Since this temperature is much lower than 
the spin gap, all the signals are assigned to be the exci- 
tations from the ground state. In this section, we show 
important features of an one-triplet excitation:a creation 
of single triplet in a dimer, and discuss the ESR se- 
lection rule considering Dzyaloshinsky-Moriya(DM) in- 
teractions and a staggered field. To avoid confusion, 
we attach subscripts :GS (ground state) and BS(bound 
state or multiple-triplet excitation) to the terms, sin- 
glets and triplets. For examples, singletcs, singletss 
and tripletss represents the singlet ground state, the 
singlet excited bound state and the triplet bound state, 
respectively. 

3. 1 Experimental features 

Figure 2 (a) shows a typical temperature dependence 
of ESR spectrum. The intensity of a broad peak at 13.3 
T decreases rapidly below 10 K. It indicates that the 
peak is caused by the transitions of the thermally ex- 
cited triplets as shown by dotted arrows in the inset. At 
1.6 K, two weak signals marked by arrows are observed. 
Since the intensity of the weak signals increases with 
decreasing temperature, these are assigned as the transi- 
tions between the ground singlet and the excited states, 
which is indicated by the solid arrows in the inset. While 
the similar temperature dependence is found for B||a in 
Fig. 2 (b), differences are found in the splitting and the 
intensity of the two peaks. As is well known, an ESR 
transition between the ground singlet and the excited 
triplet states due to the mechanism of magnetic dipo- 
lar transition is usually forbidden.^''^ A non-secular term 
such as DM-interaction, as proposed recently by Cepas 
et al, may be the leading origin of the break down of the 
selection rule.^®^ This point will be discussed later. 

In Figs. 3 (a) and 3 (b), ESR spectra measured at 
different frequencies are collected together for B||c and 
B||a, respectively. In these figures, each ESR spectrum 
is plotted in such way that its background coincides with 
the measured frequency, so that the resultant plot can be 
regarded as a kind of the frequency-field diagram which 
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Fig. 2. Temperature dependence of ESR spectra at 428.6 GHz (a) 
for B\\c and (b) for B\\a. The sharp peak at 15.3 T is the signal 
of DPPH used as a field maker. The inset shows a schematic 
energy diagram of triplet with a spin gap. 



Fig. 3. Summarized magnetic excitation spectra spanning 0-40 T 
and 0-1100 GHz at 1.6 K, (a) for B\\c, (b) for B\\a. To reduce 
the data volume, the data points are averaged with the interval 
of 0.1 T. 



furthermore show intensity variation of excitations. A 
detailed frequency field diagram is obtained by tracing 
different types of peaks in the plot. Frequency- field dia- 
grams thus obtained are depicted in Figs. 4 (a) and 4 (b) 
for B\\c and for B\\a, respectively. A contour map of the 
intensity in the frequency- field plane as shown in Fig. 5. 
Those figures are useful to see the overview of the excited 
statc^ cnc^rgy spectra in field-frequency plane. A compli- 
cated structure of the spectrum is noticed, that cannot 
be expressed well by a simple frequency-field plot. In the 
following, we show important features of the one-triplet 
excitation using the set of these figures. 

In Figs. 4 (a) and 4 (b) , it is noticed that two low- lying 
triplet excitations:Oi and O2 exist. The values of the 
zero field gaps are close to the spin gap estimated from 
magnetic susceptibility and high field magnetization'^' 
and thus these modes are assigned as one-triplet excita- 
tion. More precisely, the spin gap Eg is determined to 
be 722±2 GHz, which is evaluated as the average of the 
zero field intercepts of Oi and O2. The splitting between 
Oi and O2 is largest for B\\c and the zero field energy 
gaps of Oi and O2 are 764±2 GHz and 679±2 GHz, 
respectively. The existence of two modes is generally ex- 
pected when a unit cell contains two magnetic ions. The 
two modes, originally degenerated, may be split by an 
anisotropic inter-dimer interaction and/or other pertur- 
bations. In the present compound, the splitting of the 
two-modes is interpreted fairly well by including a inter- 
dimer DM-interaction.^^) 

When a magnetic field is rotated from the c-axis to 
the a-axis, the splitting shows a monotonic decrease and 
takes a non-zero minimum for B\\a as we reported in 
the previous work(see Fig. 4 of Ref.[15]). The angular 
dependence of the upper-going branches of Oi and O2 
is more clearly shown in Fig. 6. The two modes are 
almost parallel to each other for B\\c, while they be- 
come very close each other for B\\a accompanied by the 
bending around zero field. These behaviors indicate that 
the principle axis of the interaction causing the splitting 
points to the c-axis. 

It should be stressed here that a small but intrinsic 
splitting exists between Oi and O2 for B\\a. It should be 
noted that the splitting is not caused by the sample mis- 
orientation or anisotropy of 5-value. This can be partly 
confirmed if one measures the two modes for B||[110]. 
In this case, as shown in Fig. 1, the magnetic field is 
parallel to the A-dimers and is perpendicular to the B- 
dimers. A larger splitting would be expected, if it was 
caused by the difference of g- value between the parallel 
and the perpendicular directions of a dimer. However, it 
is found that the splitting is nearly identical with that 
for B||a(not shown here). The result indicates that the 
splitting observed for B\\a is not caused by the g- value 
anisotropy. 

The essential difference between Oi and O2 in B\\a 
is found in the following mixing behaviors. Firstly, it is 
found that the magnitude of the mixing with singlctss 
states is different between Oi and O2. In Fig. 4 (b), Oi 
shows a break around 860 GHz, while O2 shows another 
break around 650 GHz. These breaks may be caused by 
an anti-level-crossing with an excited singlet state(this 
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Fig. 6. Angular dependence of Oi and O2 in the oc-plane. 9 
is the angle between the a-axis and an external field. Lines are 

eye-guides. Anisotropy of g-valuc is plotted in the inset, which is 
measured at 25 K with 135 GHz. The negative angle corresponds 
to the rotation of an external field from the a-axis to the 6-axis. 



Fig. 4. Frequency field diagram at 1.6 K (a)for B\\c and (b) for 
S||a. O, T, Q, B, M are one-triplet, tripletg^, quintetss, 
new signal and paramagnetic like signal, respectively. For de- 
tails, see text. The subscript W denotes weaJj signals which are 
not well classified. Solid line shows the paramagnetic resonance 
line(j^ = gfiBH, i/:frequency, g:g-valuc). Dashed lines are eye 
guides. Magnetization curves taJsen from Ref.[20] are plotted 
together. 



Fig. 5. Contour plot of ESR absorption intensity for B\\a at 1.6 
K. The color scale is arbitrary. 



point will be discussed later). Secondly, it is noticed that 
Oi and O2 round very differently around He where a mix- 
ing with the ground state occurs. The mixing is strongly 
related to the symmetry of corresponding modes, and 
thus these findings indicate that the symmetry of Oi is 
different from that of O2. 

Next we discuss the intensity of the one-triplet exci- 
tation. For B\\c, a clear difference is found in the field 
dependence of intensity between Oi and O2 as shown in 
Fig. 3 (a). The intensity is nearly field independent for 
Oi, while it is sensitive to the magnetic field intensity 
for O2. Namely, the intensity of the down-going branch 
of O2 increases with increasing field and that of the up- 
going branch shows an opposite field dependence. In 
addition to this, a rapid decrease of intensity is observed 
for both branches around zero field. The behavior sug- 
gests that the origin of the break down of the selection 
rule is different between Oi and O2 in B\\c. For B\\a, 
the behavior of Oi and O2 is much different as follows. 
First, the intensity is nearly field independent for both 
Oi and O2. Secondly, the intensity is stronger than that 
in B\\c. The second feature suggests that a non-secular 
term causing these transitions should be larger for B||a. 

We compare the present result with that of the FIR 
spectroscopy.^^) It should be reminded that the defini- 
tions of absorption intensity are not exactly same each 
other. In our case, the intensity is normalized by the 
intensity of DPPH used as maker or by the intensity of 
paramagnetic signal. In the FIR experiment, the mag- 
netic intensity is obtained by taking a difference between 
a high temperature spectrum and a low temperature one. 
In spite of this difference, two results are roughly consis- 
tent each other. In Fig. 4 of Ref. [12], a field is applied 
along the a-axis. The intensity decreases with increasing 
field for the c-polarization, while it increases for the ab- 
polarization. In case of ESR, un-polarized radiation is 
used and thus the intensity should be averaged over the 
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two polarizations. If we average the data of different po- 
larizations obtained by FIR spectroscopy, it can happen 
that the intensity depends on magnetic field only weakly 
for the cancellation, which agrees with the ESR result 
for B\\a. 

3.2 Selection rule and DM-interaction 

In the following, we discuss possible interactions al- 
lowing the finite intensity of the onc-triplct excitations, 
which is originally forbidden. We consider an intra-dimer 
DM-interaction: Dd, the c-axis and the afe-plane compo- 
nents of inter-dimcr DM-interaction: {Dc and Dab, re- 
spectively) and a staggered field, all of which are allowed 
to exist by its crystal structure. Wc do not consider 
a novel mechanism:dynamical DM-interaction proposed 
recently.^®' ^^-^ Now the point to be clarified is what as- 
pect of the experimental result can be interpreted with 
those four terms. We would refuse to deny or confirm the 
existence of the dynamical DM-interaction because the 
judgement cannot be made with present experimental 
results. 

We employ the selection rule established by Sakai, 
Cepas and Ziman for DM-interaction. Since the 
present experiment is made in Faraday configuration, one 
simple rule is applicable as follows. 

DM-vule: Field independent intensity is observed only 
when the DM-interaction has a component along the ex- 
ternal magnetic field and the intensity is quadratic to the 
component. 

Among three possible DM-interactions, Dc exists irre- 
spective of the small buckling of the edge-shared CUO4 
planes as shown in Fig. 1. Dj, and Dab exist only when 
this buckling is taken into account. At the same time, 
a staggered field bcc;omcs possible for the alternation of 
a principle- axis of g-tensor as shown in the inset of Fig. 
1. Considering the large anisotropy of (/-value in the 
ac-plane(ga=2.05 and gc=2.28) and the small one in the 
ab-pla,ne{Ag < 0.01, see Fig. 6 inset), a sizable staggered 
field is expected only for J5||c. Those three terms caused 
by the buckling may be small, bcc;ause of the small tilt- 
ing angle. An experimental estimate of Dab is made 
by Kakurai et al.^^^ and it is found that the magnitude 
of Dab is about 30 % of that of D, - This estimation, 
however, looks overestimated since the buckling angle is 
quite small. The important point is that the practical 
effect is observed with such small buckling and thus Dc 
cannot be neglected. It should be also noted that Dd 
is much effective than Dab and Dc- It is simply be- 
cause D(i operates directly on the two spins in a dimcr. 
From the crystal symmetry it is evident that D^ lays 
on the a6-plane and is perpendicular to the direction of 
a dimer bond. Since two kinds of dimers are orthogo- 
nally arranged in the a6-plane, the sum of the squares of 
a field-parallel-component of D^ over the A-dimers and 
the B-dimers is constant for any field direction in the 
plane. 

Let us examine, if the DM-rule is compatible with the 
experimental results. Considering the DM-rule, Dc gives 
rise to constant ESR intensity for B\\c. This expecta- 
tion is consistent with the behavior of Oi in which the 
intensity is field- independent. A precise estimation of 



Dc=0.18 meV is made from the dispersion measured by 
neutron scattering. This value is a standard when we 
estimate the magnitude of DM-interaction from the ESR 
intensities. 

Contrary to Oi, the field dependence of the intensity 
of O2 is not compatible with the DM-rule and thus DM- 
interaction is not the origin of the transition O2. As 
shown in the previous section, the intensity of O2 shows 
the characteristic field dependence. It is noticed that 
the behavior is similar to the case of a staggered field 
with Voight configuration. In fact, a staggered field 
is expected for this field orientation as mentioned above. 
However, it cannot be applicable to the present exper- 
iment because the Faraday configuration is used. It is 
clear that neither DM interaction nor a staggered field 
can cause the transition O2 in the lowest order from the 
selection rules. A possibility is that a combined higher 
order term may cause a finite ESR intensity. The fact 
that the intensity of O2 is much weaker than that of Oi 
is consistent with the contribution of the higher order 
term. A similar situation is found in CuGeOs, in which 
both optical and acoustic branches are observed in ESR. 
This behavior cannot be interpreted by the lowest or- 
der term of DM-interaction. ^'^^ As discussed above, it is 
rather plausible that a combined higher order term of 
DM interaction and a staggered field gives the transition 
O2. A theoretical calculation of the intensity is necessary 
to clarify the present proposal. 

Next we discuss the case of B\\a. Since both Dj, and 
Dab have the field parallel components, the field indepen- 
dent intensity observed in the experiment is compatible 
with the DM-rule. The fact the intensity is larger in B||a 
than that in B\\c indicates that the field parallel DM- 
component is larger in B\\a. As mentioned above, the 
square of the magnitude of Dab is one-order smaller than 
that of Dc and thus Dab causes the minor contribution 
to the intensity. Hence, Dd should be dominant in B\\a. 
Another piece of the evidence of the major contribution 
of Dd is that the intensity of Oi and O2 is almost con- 
stant when we rotate an external magnetic field within 
the afe-plane(not shown here). It is consistent with the 
fact that the sum of the squares of the field parallel com- 
ponents of Dd in A and B-dimers is constant under this 
field rotation. Moreover, D^. operates directly on each 
dimer, while Dab does not. All these facts indicate that 
the Dd is the main cause of the selection rule breaking 
in B\\a. 

Finally, we would like to make a comment on an ex- 
change anisotropy. In our previous work, we did not con- 
sider the inter-dimer DM-interaction and mentioned that 
exchange anisotropy may cause the splitting between Oi 
and O2 and the forbidden transitions. Recent proposal 
by Cepas et al. that DM-interaction is the leading term 
is more plausible than the exchange anisotropy. Now we 
consider that a simple exchange anisotropy term does 
not cause the mixing between the singlet ground and the 
excited triplet states. However, a consideration of the 
problem from more general view point may be interest- 
ing.2.32,33) High resolution of ESR can resolve even a 
small change of triplet excitation caused by these terms. 



Hiroyuki NojiRi, Hiroshi Kageyama, Yutaka Ueda, Mitsuhiro Motokawa 



6 

§4. Multiple-triplet bound states 

As shown in Figs. 3 (a) and 3 (b), many high energy 
excitations, which have been interpreted as multiple- 
triplet bound states/^^ account for a considerable weight 
in the excited state. In the present system, a sizable 
inter-dimer interaction exists. This interaction couples 
two (or three) dimers and causes a bound sate with to- 
tal spin S'=2(or 5=3). A 5=2 bound state, for example, 
consists of three sub-states; singlet, triplet and quintet. 
Usually such bound state appears as a continuum. In an 
orthogonal-dimer system, however, discrete states are ex- 
pected for the localization of triplets. In this section, we 
describe the behaviors of those bound state excitations 
in detail. 

4-1 Singlet bound state 

A singlet bs is an excited singlet state consisting of two 
interacting triplets on nearby dimers and the energy is 
given by subtracting the binding energy from the twice of 
Eg . For a frustrated spin gap system, it is expected the- 
oretically that excited singlets can dominate the states 
below Eg. In highly frustrated system, an excited singlet 
lies in the low energy range close to the ground state. In 
another word, the magnitude of the binding energy is 
a yardstick for the degree of frustration. In the present 
case, such a low energy excited singlet indicates the close 
distance from a quantum critical point. 

In the present work, an anti-level-crossing between a 
singletss and a one-triplet excitation is found as shown 
in Figs. 7 (a) and 7 (b). In Fig. 7 (a), as the fre- 
quency increases toward 646 GHz, O2 deviates from Oi 
and shifts to the low-field side with a considerable broad- 
ening. Above this frequency, O2 approaches to Oi from 
the high-field side. For Oi, a conventional successive Zee- 
man shift is found. The bending of Oi around zero field 
is caused by the anisotropic inter-dimer DM-interaction 
as discussed in the previous section. The break of O2 
shown in Figs. 7 (a) and 3 (b) is a typical behavior of 
an anti-level-crossing. Although the singlet bs cannot be 
observed directly by ESR, the anti-level-crossing exhibits 
clearly the existence of the singlet state TSi at 646 GHz. 

A similar anti-level-crossing is found at 860 GHz be- 
tween the singlet TS2 and Oi as shown in Fig. 7 (b). 
When the frequency is very close to TS2, Oi disappears 
and only a single peak of O2 is found. The magnitude of 
the anti-level-crossing at TS2 is larger than that at TSi. 
It indicates that the coupling between the singletsg and 
the one-triplet state is larger in TS2. The present result 
also exhibits that the symmetrys of the one-triplet exci- 
tations Oi and O2 arc different each other as mentioned 
in the previous section. In Fig. 7 (c), a frequency- field 
diagram is plotted, where a magnetic field is in the aft- 
plane but not parallel to the a- nor the c-axes. An anti- 
level-crossing is found in both Oi and O2 at both TSi 
and TS2(four breaks in total). This behavior shows that 
two-modes Oi and O2 arc mixed each other when an ex- 
ternal field is tilted from the principle-axes(a or c) of the 
crystal. It implies that the anti-level-crossing is closely 
related to the anisotropic part of exchange coupling such 
as DM-interaction. 



Fig. 7. (a)ESR spectra around the anti-level-crossing between O2 
and TSi (dotted line). Dashed lines are eye-guides, (b) Same for 
Oi and TS2 (dotted line), (c) Frequency-field diagram at 1.6 
K, where an applied magnetic field is tilted 15 degree from the 
a-axis in the ac-plane. Dashed lines are eye-guides. 



Let us compare the present results with the Ra- 
man scattering experiments, where several singlet bound 
states are found. A clear peak was observed at 30 
cm~-'^(900 GHz) in Bi polarization(sec Fig. 1 of the 
Ref.[8]). In B2 polarization, a peak was also found at 30 
cm-i with a shoulder at 24 cm" ^ (720 GHz). The shoul- 
dc;r appeared also in Bi polarization and a small low 
energy shift was found at 6 T. Consequently, the shoul- 
der was assigned as one-triplet excitation being weakly 
allowed for some defects. The energy of TS2 is very close 
to 30 cm~^ and thus TS2 can attribute to the 30 cm~^ 
Raman peak. Although the energy of TSi is not far 
from the 24 cm^^ shoulder found in the Raman scatter- 
ing, these two modes are different each other because the 
former is clearly the singlet and the latter was assigned 
as one-triplet excitation. Consequently the observation 
of TSi reveals that the lowest singlet bound state is lo- 
cated at 646 GHz( 21.5 cm~^). It is notable that the 
energy is lower than the spin gap of the one-triplet exci- 
tation. This fact indicates that the system is rather close 
to the quantum critical point where the energy falls to 
zero. 

Let us compare the present results with the theoret- 
ical investigations on various bound states. ^^'^^^ As 
discussed by Totsuka, Miyahara and Ueda, the represen- 
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tations E and B2 are infrared(also ESR) active and those 

A2 and Bi are infrared inactive in D2d symmetry. The 
one-triplet belongs to the £^ representation. In Ref. [35], 
the lowest singlets^ is found to have the A2^Bi symme- 
try (additional small splitting is found in Ref. [14]). The 
second lowest singletsg coming from the third neighbor 
pair in Ref. [35] has an infrared active Ai^B2 symmetry. 
Let us assign tentatively the lowest energy mode and the 
second lowest energy mode to TSi and to TS2, respec- 
tively. It is natural to assume that the ono-triplct shows 
the stronger mixing with the infrared active mode. The 
strong/weak anti-lcvcl-crossings of TS2/TS1 are consis- 
tent with the fact that the lower singletB5(TSi) is in- 
frared inactive and that the upper one(TS2) is infrared 
active. However, since only two modes are found in the 
present experiments, the definite assignment is difficult. 
A theoretical investigation by considering the mixing ef- 
fect is needed for further understanding. Such a calcula- 
tion may be useful to examine the wave function of the 
bound states. 

4-2 Triplet bound state 

In Fig. 8, an example of ESR spectrum is shown for 
i?||a. Besides onc-triplct excitations Oi and O2, eight 
peaks are found. These additional peaks are easily clas- 
sified into triplctsBs(Ti^T5) and quintctSBs(Qi~Q3) 
by Zeeman shifts in the frequency-field diagram. In the 
inset of Fig. 8, temperature dependence of intensity is 
plotted for different peaks. Below 7 K, the intensity in- 
creases with decreasing temperature with similar tem- 
perature dependence for all three peaks. The behavior 
confirms that the additional peaks are the excitations 
from the ground state. 



Fig. 8. Anisotropy of ESR spectrum at 1024 GHz and at 1.6 K. 
Dotted arrows indicate the splitting of Oi and 02- Lines connect 
corresponding peaks. Temperature dependence of intensity is 
shown in the inset. See text for the notation of peaks. 



A most significant feature in Fig. 8 is the sizable in- 
tensity of the bound state excitations. It is the evidence 

of the large suppression of the one-triplet process in the 
present system. It also supports an idea that a bound 



state, propagating with correlated motion, has consid- 
erable spectral weight. ■^^^ In general, ESR intensity is 
not a spectral weight itself and is mediated by transition 
matrix elements mixing the ground and excited states. 
However, the intensity can be compared directly between 
the one-triplet state and the triplet's . The argument is 
based on the fact that a total spin number of the system 
changes by one in both types of excitations. From this 
basic view point, we notice that the magnitude of matrix 
element may not be so different (at least in the same or- 
der) between these two processes. In fact, the intensity 
ratio of O1+O2 to Ti is comparable with that of the first 
excited peak to the second one in neutron scattering (see 
Fig. 4 of Ref.[16]).^^' This fact supports our presump- 
tion that ESR intensity of the triplet states(in both the 
one-triplet and the bound state) roughly corresponds to 
the spectral weight. 

Another important finding is that the bound state 
spectrum splits into many sharp discrete peaks instead 
of a conventional two- triplet continuum. It suggests that 
many different types of bound state is formed for dif- 
ferent interactions acting between the distant triplets on 
dimers. The fact that these different states do not merge 
into a single band suggests that the dispersion of the 
bound states is weak. 

Let us evaluate the zero field energy gap of each bound 
state. We simply evaluate a gap by extrapolating each 
mode with a linear equation. This procedure may cause 
non-negligible error, if the level of a bound state has 
a large bending around zero field as is observed in Oi 
and O2. To check this point, we compare a zero field 
gap obtained for B\\a with that for B||c and found that 
the difference is at most 10 GHz. It is much smaller 
than that of Oi and O2(~40 GHz) and is comparable 
with experimental errors. Hence, the above mentioned 
procedure is justified. The energy gaps of the bound 
states thus evaluated are listed in Table I. Wc speculate 
that the anisotropic zero field splittings of the bound 
states are reduced effectively by the hopping effect or 
by the spread of the wave function over distant dimers. 
Another finding is that the anisotropy of the intensity 
of the tripletsss except Ti is similar to that of the one- 
triplet excitations. It shows that the origin of the break 
down of the selection rule may be the same for both 
cases. 

Table I. Zero field energy gap of the bound states. The notation 
is identical to that in Figs. 4 (a) and 4(b). 



Triplet 


Gap (GHz) 


Quintet 


Gap (GHz) 


Ti 


1140 ± 5 


Qi 


1390 ± 15 


T2 


1170 ± 5 


Q2 


1560 ± 15 


T3 


1190 ± 5 


Q3 


1600 ± 15 


T4 


1225 ± 5 


Q4 


1652 ± 15 


Ts 


1350 ± 5 


Qs 


1729 ± 15 


Te 


1220 ± 5 






T7 


1264 ± 5 






Tg 


1303 ± 5 







Several ESR spectra are shown in Figs. 9 (a) and 9 
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(b) around the one-triplet gap. The spectra change dras- 
tically in a tiny range of frequency with both the shifts 
and the intensity changes. This behavior may be caused 
by the increasing mixing between the tripletss and the 
one-triplet state. Most spectacular feature is found at 
736.1 GHz for both B\\a and B\\c. A broad and in- 
tense peak, including many different sub-peaks is found 
in high field side. Since the anomaly is found only at the 
identical frequency irrespective of field direction, it may 
be caused by the mixing with a horizontal level lying at 
this particular energy. Considering the no-mixing of the 
triplet_BS with TS2(no break is found in Figs. 4 (a) and 
4 (b)) and the no-mixing of the one-triplet states with 
the new horizontal level, it may not be a singlet bound 
state. The 5z = branch of the one-triplet of zero wave 
vector point is located at 722 GHz and it cannot be the 
candidate of new singlet state. It is found that, at the 
(1.5, 0, 0) point, an ^2 = branch of the one-triplet is 
found at exactly at this energy by neutron scattering(see 
Fig. 3 of Ref.[30]). This Sz =Q branch may be responsi- 
ble for the anomaly at 736.1 GHz. Such mixing between 
the excitation at zero-wave vector points ( (0, 0, 0) and 
equivalent points) and that at 7r-wave vector point such 
as (1.5, 0, 0) is possible by some non-secular terms such 
as a staggered field. Weaker but similar anomalies are 
found in Fig. 5 as horizontal high intensity lines, which 
may be caused by the similar origin. 



Fig. 9. ESR spectra around the one-triplet gap at 1.6 K (a) for 
B\\c and (b) for B\\a. Dashed and sohd hnes are guides of Ti 
and T4, respectively. Vertical bars show Oi and O2. B is a peaJs 
appearing around H^- 

Let us compare the present results with theoretical 



calculations. As shown in Table I, all eight bound 

sates are below the bottom of the two-triplet continuum 
(2Eg—14:4:4: GHz) aud thus these states are stable as a 
two-triplet bound state. In the reference p^-*], only five 
modes (reducing to three with degeneracy) are found be- 
low the bottom. All the mode are infrared active and 
thus we can attribute at least the five from the all eight 
modes to the theoretical ones. In the exact diagonal- 
ization result which is more reliable around a quantum 
critical point, two peaks are found side by side around 5 
mcV(1208 GHz). The energy is close to those of Ti and 
T2. In other theoretical works, the number of states is 
also less than that fotmd in the present experiment.^"' 
To have more bound states, the interactions between dis- 
tant dimers should be taken into account (an inter-layer 
coupling may also contribute). Finally we mention the 
instability of a triplet bs due to the transition to another 
ground state, proposed in theoretical work.^°) In spite 
of detailed experiments, no such anomalous behavior is 
found in the present work. 

4-3 Quintet bound state 

In the following, we show the first observation of a 
quintet bs in the system. The nature of the magnetiza- 
tion appearing around the He is also discussed. As shown 
in Fig. 8, the intensity of a quintetss is weaker roughly 
one order of magnitude, which is due to the higher order 
excitation process of the total spin change:AS'=2. Five 
quintetss modes are found in the present work, which 
are easily identified for the typical double Zeeman shift. 
A zero energy gap of a quintet bs has been extrapolated 
from the experimental data. Since the number of data 
points is limited, an error of the gap is about 15 GHz as 
listed in Table I. The anisotropy of the gap is not found 
within the experimental error. 

It is notable that only Qi is below the two-triplet 
continuum threshold starting at 2_Eg=1444 GHz, while 
the other four modes are located above the threshold. 
Since the energy of Qi is below the threshold, it is 
assigned as a quintetss made up of two-triplet parti- 
cles (four spins). Let us discuss the nature of other high 
energy quintetssg. In the isolated dimer limit of four 
spins, the tripletsg and quintetss are separated by the 
eflfective antiferromagnetic interaction acting among the 
four spins. In this limit, the higher energy of Q2 ~Q5 
can attribute to this interaction and those peaks are con- 
sidered as two-triplet bound state. This interpretation 
has a difficulty because a multi-particle bound state is 
unstable when the energy is above the corresponding in- 
coherent continuum in general. It is totally unclear if the 
triplet and the quintet sectors of the bound state can sep- 
arate clearly from the continuum spreading above 2£Jg. 
A more natural idea is to assign Q2 ~Q5 to three-triplet 
bound state and in fact such three triplet state is found 
in theoretical calculations, which is shown in later. The 
existence of three-triplet bound state shows the signifi- 
cant contribution of higher order excitations. 

The existence of Qi below the threshold causes the 
peculiar feature that the quintetss touches the ground 
state before the closing of one-triplet gap. In spite that 
a finite mixing between the ground and excited states 
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modifies this simple situation to the more complex one, 

a characteristic contribution of a quintet is found exper- 
imentally. The enlarged spectra of Figs. 4 (a) and 4 (b) 
around the He are shown in Figs. 10 (a) and 10 (b). 
It is noticed that ESR modes bend around the and 
the behavior is much different between B\\c and B\\a. 
In fact, the antistory of Qi cannot be normalized by the 
g-value and thus is originated by the anisotropic mixing 
term. This anisotropy of ESR spectra should correspond 
to the anisotropy of the magnetization curves. It is be- 
cause the magnetization around the He is caused by the 
mixing between the ground and the excited states. In 
another word, the magnetization should be dominated 
by the wave function component of the modes enhanced 
by the mixing. Keeping this point in mind, let us exam- 
ine the magnetization and the frequency-field diagram 
carefully. 

For B\\c, the magnetization becomes finite above 15 
T and at the same time O2 starts to deviate from a 
straight line. This fact shows that the one-triplet com- 
ponent dominates the magnetization. The intensity of 
Qi is enhanced around 20 T presumably for the weak 
anti- level-crossing with Oi. Around 21 T, Qi disappears 
and this fact indicates that the quintet component in the 
magnetization is small. After the anti-level-crossing, Oi 
bends considerably and loses the intensity gradually. It is 
noticed that O2 shows a turn around like behavior above 
20 T, which indicates the strong mixing with the ground 
state. On the other hand Oi goes down again above 26 
T and disappears suddenly at the jump of magnetization 
to the g-plateau as shown in Fig. 4 (a). The origin of the 
flat part of Oi around 25 T might be a anti-level-crossing 
with a new singletss, however, a further investigation is 
necessary for this point. The strong mixing of O2 with 
the ground state indicates the dominant contribution of 
the one-triplet component to the magnetization. It is 
also noticed that Q3 is cinhanced with a large round- 
ing above 20 T, which shows the mixing of three-triplet 
state with the ground state. The contribution of differ- 
ent components to the magnetization indicates that the 
magnetization is not uniform and homogeneous. It is 
much different from a usual spin gap system in which 
a magnetization consists of homogeneously distributed 
triplets. This characteristic behavior may be caused by 
the localization of excited states in the present system. 
Recently, a homogeneous internal field is found by NMR 
below the g-plateau, however the difference is presum- 
ably due to the low observation frequency in NMR, in 
which the distribution of an internal field is averaged.''*^ 

Let us discuss the case of B\\a shown in Fig. 10 (b). 
While the behavior of Oi is similar to that in B||c, both 
Qi and O2 continuously go down toward the ground 
state. The enhancement of Qi and the disappearance 
of O2 indicates that the quintetss component dominates 
the magnetization. Other higher energy quintetss states 
also exhibit enhancements and bendings, which indicate 
the strong mixing of the quintetsss with the ground 
state. An observation of a small break of O2 in i3||a(see 
Fig. 4(b)) and at 22 T may be caused by the anti-level- 
crossing with unknown singletss at 100 GHz, which is 
also found in Bile. 



Fig. 10. The low energy part of ESR spectra around He (a) for 
B\\c and (b) for B\\a. The dashed lines are eye-guides. The 
vertical bar denotes Qi. 



As shown above, the magnetization below the g- 
platcau contains different types of stat(^s and does not 
consist of homogeneously distributed triplets. A domi- 
nant contribution of a quintetsg is found for B\\a. It 
is an interesting point, if a quintetsg coupled by an ef- 
fective repulsive interaction shows a condensation at H^, 
which was proposed by Momoi and Totuska.^^^ Experi- 
mentally the domination of the quintetss is found clearly 
in the present work. However, it may be a quenched and 
localized state rather than a condensation. In the case 
of a condensation, a magnetization should not contains 
a several different states as observed in the present case. 
The localization of excited states and the resulting heavy 
effective mass of a quintet 55 may prevent the condensa- 
tion of a quintetss. 

Finally we compare the present results with the the- 
oretical calculation of the quintetsg. In the reference 
P^^], the lowest mode has the A2®Bi symmetry, which 
is infrared inactive. Since a quintetss shows very weak 
intensity, an infrared inactive mode should not be ex- 
cluded from the candidate. The fact that the number of 
two-triplet quintetss is only a few is consistent with the 
present result. However, a large number of a three-triplet 
quintetss found in the experiment calls for further the- 
oretical investigations. 

§5. High Field Phase 

In this section, we discuss ESR signals observed above 
He and their concerns with the magnetization plateaux. 
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In this region of magnetic field, the most of the inten- 
sity concentrates into the signals B and M as shown in 
Figs. 3 (a) and 3 (b). These two signals alternate at 
the magnetization jump to the ^-plateau as shown in 
Fig. 11. The position of M is identical with that of 
the corresponding paramagnetic resonance and the ex- 
trapolation of M crosses the origin of the frequency-field 
diagram for both B\\c and B\\a as shown in Figs. 4 (a) 
and 4(b). On the other hand, B shows the asymptotic 
behavior to the paramagnetic line. When a spin gap is 
closed, a 3D ordercd-state or a disordered ground state 
dominated by gap- less triplets appears. In the former, a 
field variation of a spin wave is observed as an antiferro- 
magnctic resonancc(AFMR). In the latter, the Zecman 
splitting of a triplet causes ESR at an identical reso- 
nance field with that of a paramagnetic resonance. In the 
disordered state, the strong correlation among the spins 
does not cause a shift of resonance field. It is because 
the spin operator commutes with the isotropic exchange 
coupling. "^^^ 

The behavior of M shows that the signal is the ESR 
of a disordered state. It is noticed that M continues at 
j-plateau and a low-energy excitation is absent. Such 
low energy discrete excitation is observed at ^-plateau, 
which corresponds to the creation of a new triplet into 
a super-lattice structure at a plateau. These features 
can be understood as the analogy of Ising spin cluster 
excitation. When Ising anisotropy is added gradually to 
a Heisenberg spin system, a spin wave excitation shows 
a crossover to a discrete spin cluster excitation. In the 
present system, Ising anisotropy corresponds to the hard- 
ness of a plateau. If the plateau is soft, a creation of 
single triplet (spin cluster like excitation) does not occur 
and a perturbation by the radiation transmits through 
the transverse components of spins as a wave. This ex- 
citation is similar to a spin wave excitation because it is 
likely a gap-less Goldstone mode. Momoi and Totuska 
proposed that a "supersolid" state appears between the 
plateaux, in which a spin density wave and superfluid 
of spins coexist. The latter component gives rise to 
a gap-less ESR. The present result indicates that a su- 
perfluid component (gap-less mode) is significant at \- 
plateau and that the plateau is not completely static at 
1.6 K. Recently, the static super lattice structure is found 
by NMR in the g-plateau at 30 mK.^*) For i-plateau, a 
static super structure may also appear if we measure in 
much lower temperature below 1.6 K. This point remains 
for further experimental investigations. 

Although the behavior of M is rather conventional, 
B shows unusual features:(l)location above the param- 
agnetic line, (2)non-linear field dependence, and (3)dis- 
appearance at the jump to the g-plateau. Consider- 
ing these features, we propose two possible origins of 
B:(a)AFMR, (b)transition 5^ = 1 ^ 5^ = in the low- 
est triplets. First we discuss the possibility of AFMR. 
The asymptotic behavior of B matches with that of 
AFMR mode when a magnetic field is perpendicular to 
the plane of spin laying. In AFMR, a deviation from the 
paramagnetic line is caused by anisotropy and thus the 
deviation should change for different field orientation. 
However, the deviation of B is nearly equal between B\\c 



and B||a. The most fundamental diflaculty of the inter- 
pretation by AFMR is that no evidence of 3D-ordering is 
found below the ^-plateau. Actually, in NMR measure- 
ment, no static internal field is found in this field range, 
which shows non-existence of 3D-order.'^®^ A possible in- 
terpretation is that a quasi static order takes place and it 
fiuctuates with the frequency between ESR and NMR. In 
another word, the system is considered to be dynamical 
in the time window of NMR and to be static in that of 
ESR. However, it is not clear if such situation is realized. 

Next we discuss the possibility of the case (b). In 
the most simple picture of a isolated dimer system, the 
5*2 = 1 level is below the singlet ground state above 
H,. as shown in the inset of Fig. 11. The transition 
5*2 = 1 ^ 5*2 = shown by the arrow(i^2) represents ESR 
of a disordered magnetic state and the resonance field 
is identical with that of a paramagnetic resonance. In 
the present system, the level scheme is modified for the 
anti-level-crossing as shown by the dotted lines. The 
transition 5z = 1 ^ 5*2 = also changes to that of the 
dashed arrow(z^i). The break by the anti-level-crossing 
causes the deviation of the transition from that of a 
simple triplet. The signal B may be assigned as the 
transition vi. Qualitatively, the features (1) and (2) 
mentioned above are consistent with this interpretation. 
The strong intensity of B is also consistent with the pic- 
ture because it is the allowed transition from the ground 
state. However, a quantitative agreement is difficult 
with the simple anti-level-crossing model of the inset. 
In the simple model shown in the inset, the relation 
vi=V2+Vd, should hold at each field, where is the tran- 
sition Sgs = — > = l(note that arrows in the inset 
are shifted to prevent the overlapping). It is because of 
the anti-symmetry of the anti-level-crossing around the 
He- However, this relation is not found among Oi, O2 
and B as can be seen clearly in Figs. 4 (a) and 4 (b). 
Namely the deviation of B is rather independent of the 
field direction, Oi and O2 show the anisotropic behav- 
iors. As shown above, both proposal are not sufficient to 
explain the behavior of the signal B. The nature of the 
state below the ^-plateau seems to be unconventional, 
which requires further microscopic investigations by us- 
ing different types of probes. 

§6. Summary 

To be summarized, the following important features 
are found by the detailed investigation of the magnetic 
excitation spectrum up to 1.3 THz and up to 40 T. 
(l)Difference of the symmetry between the nearly de- 
generated one-triplet excitations Oi and O2 are found in 
the field dependences and in the anti-level-crossing with 
the singletss and the ground state. 

(2) Two-singlet bound states are found by the anti- 
level-crossing with the one-triplet excitation. The one 
of those is located below the spin gap of the one-triplet, 
which indicates the proximity to the quantum critical 
point of the present system. 

(3) Several triplet bound states are found and the zero 
field energy gaps are all below the two-particle contin- 
uum threshold. 
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Fig. 11. The ESR spectra around the g-plateau for B\\c at 1.6 
K. The thick and thin dashed lines indicate the paramagnetic 
line and the field of the magnetization jump to the i-plateau, 
respectively. The dotted line is an eye-guide for the signal B. The 
inset shows the schematic energy level of an isolated S = 1/2 
antiferromagnetic dimer. Sqs = is the singlet ground state. 
For others, see text. 



(4) The gap of the lowest quintetss is found to be be- 
low the two-particle continuum threshold and it crosses 
with the one-triplet excitation before the closing of the 
normal one-triplet spin gap. Some of the quintetss 
modes are assigned clS cL three-triplet bound state. 

(5) The finite magnetization appearing below the g- 
plateau contains various types of magnetic components 
caused by the anisotropic mixing and the content changes 
with the field direction. 

(6) Two models are proposed for the unusual ESR 
mode B observed below the magnetization jump to the 
|-plateau. 

(7) The domination of gap-less ESR signal M indicates 
that the ^-plateau is rather soft compared to the ^- 
plateau. 

(8) Qualitative interpretation is made for the selection 

rules in the one-triplet excitation by considering intra- 
dimer DM-interaction, inter-dimer DM-interaction and 
a staggered field. 

The present work exhibits various novel features of 
higher order multiple-particle excitations in a quantum 
spin gap system. 
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